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Abstract. Decaying dark matter cosmological models have been proposed to remedy the 
overproduction problem at small scales in the standard cold dark matter paradigm. We 
consider a decaying dark matter model in which one CDM mother particle decays into two 
daughter particles, with arbitrary masses. A complete set of Boltzmann equations of dark 
matter particles is derived which is necessary to calculate the evolutions of their energy den- 
sities and their density perturbations. By comparing the expansion history of the universe 
in this model and the free-streaming scale of daughter particles with astronomical observa- 
tional data, we give constraints on the lifetime of the mother particle, T -1 , and the mass 
ratio between the daughter and the mother particles mv/myi. From the distance to the 
last scattering surface of the cosmic microwave background, we obtain > 30 Gyr in the 
massless limit of daughter particles and, on the other hand, we obtain mj) > 0.97?tt,m in 
the limit T" 1 — > 0. The free-streaming constraint tightens the bound on the mass ratio as 
(r-7lO- 2 Gyr) < ((1 - m m /m M )/W~ 2 y 3/2 for T^ 1 < H^z = 3). 
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1 Introduction 

There are numerous astronomical phenomena which indicate the existence of cold dark matter 
(CDM), such as flat rotation curves of galaxies, anisotropies of cosmic microwave background 
(CMB), and measurements of gravitational lenses [1]. However, we know little about the 
nature of cold dark matter. In fact, at sub-galactic scales, there are some discrepancies 
between CDM model predictions by numerical simulations and astronomical observations, 
such as the abundance of substructures in dark matter halos and the inner matter distribution 
in galaxies [2, 3]. 

The discrepancies have stimulated numerous proposals to rescue the standard CDM 
model both from astrophysics and particle physics communities. One of the astronomical 
explanations for these discrepancies is gas feedback: the heated gas by astrophysical processes 
such as supernova explosions may alter the matter distribution at the inner core of galaxies 
and/or hinder the galaxy formation in sub-halos. Explanation from particle physics, on the 
other hand, includes decaying dark matter, interacting dark matter, warm dark matter and so 
on (for a review, see [4]). Among them, Cen [5] proposed the decaying dark matter model to 
solve both the overproduction problem of dwarf galaxies and the over-concentration problem 
of the inner core. 

There are many kinds of studies about decaying dark matter, in particular, constraining 
its lifetime T _1 [5-17]. For example, the light element abundances predicted by Big-Bang 
Nucleosynthesis set the limit on the lifetime if the lifetime is relatively short [12, 16]. For 
decaying dark matter with much longer lifetime, constraints can be set from the other cos- 
mological observations, such as anisotropies of CMB, the abundance of clusters of galaxies, 
halo mass profiles and so on [6, 8, 13-15]. Most of these constraints have been obtained with 
an assumption that dark matter decays into massless daughter particles. Recently, however, 
a decaying dark matter model is proposed in which the mass difference between the mother 
and daughter particles is very small and the lifetime of the mother particle is as long as 
the age of the universe [7, 18]. By using N-body simulations of galaxy halos, Peter [7] put 
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constraints on T and the mass of daughter particles through the recoil speed of daughter 
particles in the case where the mass difference of the mother and one of the daughter par- 
ticles is small. There are several theoretical proposals for models with decaying dark matter 
(see, eg., [17, 19-22]). 

In this paper we extend previous works [6, 13] by taking into account the finite mass 
effects of daughter particles with arbitrary masses. We derive a set of Boltzmann equations to 
describe this decaying dark matter scenario, and study the effects of the decay on cosmology. 
In our formulation, the masses of the mother and two daughter particles, mM,mDi, an d 
mD2i can be set to arbitrary values. For an example we assume that the mother particle is 
Next Lightest Supersymmetric Particle (NLSP), which belongs to Minimal Super symmetric 
Standard Model(MSSM). If mother is a NLSP, its mass should be larger than at least 46 GeV 
[23] . Due to the R-parity conservation, if the mother particle is a SUSY particle, one of the 
daughter particles is also a SUSY particle and the other is a standard model particle. We 
only consider the possibility that the created standard model particle is a lepton, which can 
be considered effectively massless in the decay process. In this paper, therefore, when we give 
constraints on the lifetime of the mother particle we assume that one of the daughter particles 
is massless. By solving the set of Boltzmann equations for the mother and daughter particles 
and by integrating the distribution functions we can derive the time evolution of the energy 
density and the typical free-streaming scale of daughter particles. We then give constraints 
on r _1 and the mass ratio w,di/^m from the "comoving" angular diameter distances to 
CMB and the position of Baryon Acoustic Oscillation (BAO). We also give constraints from 
the free-streaming scale of daughter particles by comparing with the Lyman a data [27, 28]. 

This paper is organized as follows. In Sec. 2, we derive a complete set of Boltzmann 
equations for mother and daughter particles to describe the evolutions of their energy densities 
and density perturbations including the decay process. In Sec. 3 we solve the background 
distribution functions of mother and daughter particles. We then show the evolution of 
their energy densities and give observational constraints from astronomical data sets on the 
lifetime of the mother particle and the mass ratio in Sec. 4. Finally Sec. 5 is devoted to our 
conclusion. Throughout this paper the speed of light c is set to unity. 

2 Formulation 

In this section, we aim to derive a set of equations for mother and daughter particles in order 
to describe the time evolutions of their energy densities and their associated linear density 
perturbations including the decay process. In order to do this, we write down a complete 
set of Boltzmann equations (both at zero-th order and first order) in this section and in 
Appendix C. The background equations in Sec. 2 will be indeed solved in our analysis to 
derive constraints from the geometry of the universe, i.e. from BAO and CMB peak positions, 
and Hubble constant today. 

To describe the decay process, we define the decay rate r((/Dj <7m)j which is the function 
describing how many daughter particles with comoving momentum grj are created for a unit 
time interval from the mother particles with momentum q^. First, let us write the Boltzmann 
equations for the distribution functions of mother and daughter particles. The Boltzmann 

1 If the condition is relaxed, it is possible that the mass of mother particle is lighter than the mass of the 
lightest particle which belongs to MSSM (see eg. ,[24-26]). However, we can still regard the mother particle 
massive enough in these cases to take the one of dauther particles as massless one. 
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equation for the distribution function of the mother particles /m(<?m) is 

dfu _ df M | dx l df M | dq M df M | drg df M _ ( dfu \ ^ ^ 

dt dt dt dx l dt dqM dt drii \ dt J c 

where rii is the unit vector in the direction of the momentum. 

The collision term in eq.(2.1) can be expressed as an integration of r(<7D> <Zm)/m with 

<?D as 

{^d¥) c = ~ I ^^fud'qv ■ (2.2) 

Similarly, the Boltzmann equations of daughter particles are 

dfpj _ df Dj | dx l dfpj | dq D df Dj | dnjdf Dj _ f df Dj \ 
dt dt dt dx % dt dqD dt drii \ dt J c 

where j is the particle index j = 1,2. 

The daughter particles in this model are created only through the decay of mother 
particles. Thus the collision term can be written as 

(^f-) c = + / T (QB,qM)fMd 3 q M ■ (2.4) 

Secondly, we derive T^d^m) explicitly. From the definition of CDM, the thermal 
motion of mother particles is so slow that we can regard them as objects at rest. To be 
specific, the kinetic energy of mother particles is much smaller than the mass deficit of the 
decay. Due to the conservation of the momentum, two daughter particles are emitted in 
the opposite directions and the amplitudes of these momenta are the same. So we can 
write the decay rate T^d^m) as a function proportional to the delta function of qo, which 
is determined by the energy momentum conservation. From the conservation, qy[ and q^ 
should satisfy the condition as 



Because clium, the solution of eq.(2.5) leads to 

1 



a2 + lM = V m Dl a2 + ^D + V m D2 a2 +9D . (2.5) 



q D ~ ^A M o 2 +B u q M ■ (2.6) 

In addition, we assume that these daughter particles are isotropically emitted. By taking 
into account the spherical symmetry of the decay, the decay rate should be written as, 

r(<?D,<7M) = ^2" <5 (on ~ \\J 1 A Md 2 + B M ql^j , (2.7) 

where T _1 is the lifetime of the mother particle and a is the scale factor. Here Ay[ and Bm 
are constants defined as 

A M = m 2 M - 2{m 2 m + m 2 D2 ) + ( m Di ~ ™ D2 ) (2 8) 

(2-9) 
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Thirdly, we consider the first order perturbations of the distribution functions /m, /d? 
(j = 1,2) as follows. It is convenient to write a distribution function as a zero-th order 
distribution, which is the background distribution, plus a perturbed function \[/, 

fu = A 0) ((?M,t)(l + ^M(^,gM,ni,t)), (2.10) 
hi = f§](QD,t)(l + * Dj (x i ,q D ,n l ,t)). (2.11) 

By substituting eq.(2.10) and (2.11) into eq.(2.1) and (2.3), respectively, and comparing the 
equations order by order, we obtain the following equations at zero-th order for the mother 
particles as 

unperturbed : =-J T(q B , gwJ/^dV = -Tf$ ■ (2-12) 

Here overdot denotes the derivative with respect to the cosmic time t. 
For daughter particles, we obtain the following equation, 

unperturbed : f§) = [ T(q D , q M )fM (QM)d 3 q M = tt- — /m^d) , (2.13) 
J J -DM <?D 



where 




(2.14) 



The first order equations are presented in Appendix C. 

The two unperturbed equations above can be interpreted as follows. The equation (2.12) 
states that all the mother particles should decay with the decay rate T. For daugher particles, 
on the other hand, the equation (2.13) means that the daugher particles with momentum 
(/D should be created at a rate proportional to r/^ (g^), where q^ is the momentum which 
the mother particles should have for the created daugher particle to have the momentum go 
after the decay. Note that because from eq. (2.13) the unperturbed distribution functions of 
two dauther particles coincide, = f$, and hereafter we simply denote them as fj^\ 



3 Calculation of background distribution functions 
3.1 Mother particle 

The temperature of CDM is low compared with the rest energy of the particle. If the mother 
particles do not decay, i.e. T = 0, the background distribution function of the particles is 
given by the Maxwell-Boltzmann function. In this case, or in the very early universe where 
the decay is negligible, the background distribution fffl is 

= (2™ m Tmo) 3 / 2 ^ (-^fk) ' (3 - 1} 

where Tmo is the present temperature of mother particles. We consider 7?im = 1.0 TeV as 
a working example. Because of the calculation performed in Appendix b, we obtain the 
temperature of mother particles as 

T M0 - 1.7 x 10" 14 K = 1.4 x 10" 18 eV , (3.2) 
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And we set the Boltzmann constant k-g = 1. In fact, the result depends only on the 



f(o) 

by the solution of eq.(2.12) as 



mass ratio mDi/mM, but not on the absolute value of myi- Including the decay, is given 



Here we have normalized the distribution function at t = as 

r+oo 

/ 47rg| I /^ ) (g M , t = 0)d<? M = 1 • 
■/ o 

The number density of the mother particles, tt-m, is determined by 

n M = Pc^dmMm > ( 3 - 5 ) 

where p c is the critical density of the universe, Qdm is the density parameter of dark matter 
normalized by the critical density at present. Note that the parameter S^dm in the above 
equation is an extrapolated value of the density parameter of mother particles without decay. 
Then the energy density of the mother particles pm is given by 



(3.4) 



Pc^DM 
PM = — 4- 



/>+oo , 

J 47r QM\/ m M a2 + ^m/m (QM,t)dq M ■ (3.6) 



m,Ma q 

3.2 Daughter particle 

In this section we derive the background distribution function of daughter particles. By 

f (o) 

'M 



substituting fj^ in eq.(3.3) into eq.(2.13), we obtain a partial differential equation for the 



unperturbed distribution function fjy as, 

/ \ 3/2 

• (0 ) _ V2T ( 1 ^ ' [Zh. ~T77_ ( ±ql-a 2 A M \ 



A" - J - ( - fess J exp ( - rt) ' m 

Since the daughter particles did not exist in the early stage of the universe, 

fSHoD,t = o) = o. 

When a mother particle decays, the amplitudes of physical momenta of two daughter particles 
are the same, and we denote it as pth- Because the physical momentum of daughter particles 
decays as oc a" 1 as the universe expands, the daughter particles which are created in the past 
should have the comoving momentum smaller thanp t h. In addition, since the thermal motions 
of mother particles are very slow compared with their mass, the time can be decided uniquely 
when the daughter particles with the momentum were created. To put it concretely, the 
redshift zj) which corresponds the redshift when the daughter particles with the present 
momentum qj) were created should satisfy 

1 1 A 

<?D = — Pth ~ z ' Am . 

1 + Z D 1 + ZD 

where Am is the mass difference between the mother and daughter particles. In the limit 
mm + m.D2 — > 0, pth is equal to Am. 
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Now let us consider the time evolution of at a fixed comoving momentum go- 
The source term of eq.(3.7) is exponentially suppressed in the very early universe when 
a 2 <C Aq^/AM, because the typical momentum of the daughter particles qn is much larger 
than the termperature of the mother particles Tmo- The source term, on the other hand, 
should be zero when a 2 > Aq^/Ayi, which comes from the energy-momentum conservation 
law. Therefore, the source is important only around t < t* where t* is defined by 

4ql-a 2 (t* q JA M =0 . 

To take advantage of the rapid convergence of the source term of /q (qn,t) we expand 
a(t) around t* D as 

<t) ~ a(t*J + a(t*J(t - i*J = a(i*J + d(t* D )e (e < 0) , (3.8) 

where e = t — t* D . Here we have omitted the higher order terms in the expansion because the 

source of (?D,i) decays exponentially backward in time for t < i* D . Then the evolution 
equation (eq.(3.7)) can be expanded as 



fS\qD,t) = fS\qD,t* qD +e) ~ o J"„ 3/2 | ) x /-a(t* qD )a(t* qD )eA M , 




and the integration of time t can be replaced with that of e. Furthermore, we can extend the 
range of integration as 

/dWO= f fg\qn,t* v +e)de~ f fg\q D , t* qD + e)de . (3.10) 

J-t* D J-oc 

This is because the term in the exponential in eq.(3.10) is very large in negative value; 
A M a(t* )a(t* ) 

( — K ) ~ — 0(10 ) <C —1 for t* around recombination, and hence the 

B M muTuo 

integration of eq.(3.9) for t = [— oo,0], i.e., e = [— oo, — t* D ] is negligible. We are then able to 
perform this integration of e analytically to obtain 



(3.11) 

We show the shape of the distribution function of daughter particles {qu^t) at present 
time in figure 1. 

We depict the distribution function of daughter particles at two different epochs, at 
matter-radiation and matter-A equalities, in figure 2. We find that the distribution func- 
tion, / D , is proportional to g£, and q D for g D < p t h/(l + z cq ) and j>th/(l + z cq ) < qr> < 
Pth/(1 + z h), respectively, where z eq and z\ are the redshifts of matter-radiation and matter- 
A equalities, respectively. These dependences can be understood as follows. First let us 
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Figure 1. Distribution function of daughter particles at present as a function of comoving momentum 
of daughter particles go, with parameters tum = 1.0 TeV, mm = 0.98 TcV, and = 0.1 Gyr. 

estimate the number density of daughter particles n. When t <C the number density of 
daughter particles is given as 

r 

^ n M (l - exp(-lY)) ~ — n M • (3.12) 
The number density tid is also expressed through an integration of as, 



n 



Jo 



Thus we can express with H and as, 

AO) n D r , . 

/d ~^~^ nM - (3 - 13) 

In the radiation-dominated epoch, the scale factor a grows proportional to i 1 / 2 and the 
Hubble parameter H is proportional to a -2 . In addition, the masses of the mother and the 
daughter particles determine the momentum of daughter particles at their creation to a fixed 
(constant) value pth- Thus qv can be written as 

qv = apth oc a . (3.14) 
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Figure 2. Snap shots of the distribution function of daughter particles (red lines) at matter-radiation 
equality t eq (left) and matter-lambda equality t\ (right). The present distribution function is also 
given (gray dotted). The green and blue dotted lines are the fitting lines with oc q^ 1 and <7 D 3 ^ 2 , which 
represent the particles created in the radiation and matter dominated epochs, respectively. We find 
that these fitting lines are in good agreement with the redlines. 



Combining the above dependences altogether, we can derive the go-dependence in figure 2 
in the radiation-dominated epoch as 

/o^D^eq) OC q^ 1 . (3.15) 

On the other hand, in the matter-dominated epoch where a oc oc a 3 / 2 , we obtain 

the distribution of daughter particles created in the matter dominated epoch in the same 
way as 

^% D ,t A )a<fo 3/2 • ( 3 - 16 ) 

We confirm that these dependencies are indeed found in figure 2. 

The energy density of daughter particles can be calculated in the same way as mother's. 
It is given by 

pd = [ +00 W D J<^ + q yS\ q ^t)d qD (3.i7) 

r+oo I 

J ^qvyjm^a 2 + q^ffi' (qr,,t)dq D ■ 



_l_Pc^DM 



m M a A Jo 



Finally the time evolution of an homogeneous and isotropic expanding universe follows 
the Friedmann equation, 



a 



a I SttG . . , . 

^— {PM + PD + PB + p<y + Pv + PA) , (3.18) 



where pb, Py, Pu and p\ are the densities of baryon, photon, neutrinos, and dark energy, 
respectively. 
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4 Result and Discussion 



4.1 Background energy densities 

By integrating the distribution functions fffl and , we obtain the time evolution of energy 
densities of mother and daughter particles, which is shown in figure 3. As is shown in the 
figure, if the dark matter particles decay, the total energy density in the universe becomes 
small comparedwith the standard A-CDM model [29, 30]. Thus the time evolution of the 
scale factor a differs from that of the A-CDM model. This leads to the different angular 
diameter distances to the last scattering surface of CMB di±{z*) and the position of BAO d z . 
The distances to CMB and the position of BAO are measured precisely by WMAP [31] and 
SDSS [32], respectively. The uncertainties of these measurements are also available in those 
papers which can be used to constrain the decay rate F and the mass ratio mm/mM, as we 
will discuss below. 




Figure 3. Evolution of energy densities of the mother (green) and daughter (blue) particles as a 
function of scale factor a with model parameters tom = 1-0 TeV, mo = 0.20 TeV, and = 0.1 Gyr. 
The total energy density of mother and daughter particles is shown as a red line. The time evolution 
for the standard CDM model is also shown for comparison (black dotted line). 

4.2 Constraints from Hubble constant, BAO and CMB 

In this paper we only consider the dark matter which decays after cosmological recombination. 
For this reason, we fix standard cosmological parameters to the values which agree with the 
A-CDM model in the early universe obtained by WMAP7 [31]. Because the decay of mother 
particles can be neglected deep in the radiation dominated era, the initial conditions of the 
dark matter energy density pyi and the scale factor a(t\) can be set as in the A-CDM model 
without decay. In this paper we use the relation which holds in the radiation dominated era, 
a(t = 0.02 sec ) = 4.60 x 10~ 12 as our initial condition. 

The "comoving" angular diameter distance to CMB, c?a (■?*), is sensitive to the deviation 
of a, where is the shift of recombination. Here d\(z*) is determined precisely by WMAP 
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project as 

d A (**) = 141161^° M pc . 

In the same sense the distance to BAO deviates from that in the A-CDM model. Ac- 
cording to the observation of BAO by SDSS [32], the BAO constraint is given through the 
variable defined as d z = r s (zj))/Dy(z), where r s {z^,) is the comoving sound horizon at the 
baryon drag epoch and zd is the redshift when photons decouple baryons. According to the 
report of WMAP7 [31], 

r s (zd) = 153.2 ±1.7 Mpc . (4.1) 
The distance Dy(z) is a function of redshift z defined as [33, 34], 

Dv{z)s (ii±4im\ ll \ (4 . 2) 



B(z) 

where d^{z) is the angular diameter distance to the point whose redshift is z. When one 
specifies a cosmological model and the evolution of the scale factor a is determined, we can 
use d z to constraint the parameters of the model. The distances do. 2 and c?o.35 are given by 
SDSS DR7 as [32], 

dtio = 0-1905 ± 0.0061 , (4.3) 
d°j» = 0.1097 ± 0.0036 . (4.4) 

We define a vector x as 

- = f " , (4.5) 

\ a 0.35 "0.35/ 

where do20'^035 are ^ ne distances based on the cosmological model to be constrained. The 
matrix C _1 , which is inverse of the covariance matrix C = (x t x\ is given by, 

1 = ( 30124 -17227\ 

" V- 17227 86977 J ' 1 ' 

where t x is transpose of the vector x. We use the value \ 2 = l xC~ l x for our chi-square test. 
Since we have two model parameters, the regions of la and 2a confidence levels correspond 
to those which satisfy Ax 2 < 2.18 and 6.30 from the minimum, respectively. 

Due to the dark matter decay, the total energy density in our model is always lower 
than that in the A-CDM model if the energy densities of dark matter in the early universe are 
fixed to the same value between the two models. Thus the Hubble parameter H of our model 
is always smaller than that in the A-CDM, which we denote as -Ha— cdm- A simple constraint 
on our model is therefore put from the current value of the Hubble parameter. The latest 
compilation determines the Hubble constant as Hq = 73.8 ± 2.4 km/s/Mpc [35] including 
systematics. By comparing the Hubble parameter with the data, we obtain constraints on 
the lifetime of the mother particle T _1 and the mass ratio mm/nriM, which are shown in 
figure 4. 

We consider the " comoving" angular diameter distance toward the last scattering surface 
(•?*)) where z* is the redshift of recombination. The distance c?a can be written as 

J H(z) J H a _cdm(z) 
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Figure 4. Constraints on the lifetime of decaying dark matter and the mass ratio between the mother 
and the daughter particles. Solid lines and dashed lines correspond to 1 a C.L. and 2a C.L. of the 
constraints, respectively. The names of observational data are also shown in the figure. 

Therefore 

d A 0*) > dA(A-CDM)0*) • (4.8) 

In figure 4 we show the constraint on the lifetime of decaying dark matter from the 
distance to BAO and the last scattering surface of CMB. The A— CDM model corresponds 
to the limits r — > and/or (1 — moi/mM) — > 0. As T becomes smaller, or moi approaches 
closer to tom, the distance cIa(z*) approaches to the value in the A-CDM model. In general, 
constraints from CMB are stronger than those from BAO. This is simply because CMB data 
are more precise than the current data of BAO. In the massless limit of daughter particles, 
we find T _1 > 30 Gyr from CMB (at la). Because we have derived the constraint from the 
background fj^\qD,t) on ly 5 this constraint is weaker than that obtained by [6, 13] in which 
the perturbations are included. At the limit T" 1 — > we obtain mDi > 0.97mM from CMB 
(at la). These values are consistent with the results obtained by the simplified estimations 
without using the distribution functions (see Appendix a). 

4.3 Free- streaming scale 

The daughter particles created by the decay of mother particles move with large velocity 
in the expanding universe. It leads that the depth of the gravitational potential of mother 
particles becoming shallower which have been created by a group of mother particles becomes 
shallower in time. In the linear theory of structure formation, the structures smaller than 
the free-streaming scale lps are erased. Since it is difficult to fully calculate the density 
perturbations of the mother and the daughter particles, here we use /fs f° r constraints on 
r _1 and uidi/ TTiM- We restrict our attention to the case where T _1 is much smaller than 
the age of universe at redshift z, T _1 <C t ~ H~ 1 (z). In this case, because almost all the 
mother particles have decayed by that time into daughter particles which are responsible for 
the depth of gravitational potential of mother particles becoming very shallow, we can place 
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Figure 5. Time evolution of the free-streaming scale of daughter particles for todi/^m = 0.98 (left) 
and 0.50 (right). The different lines in the panels correspond to the different lifetime of the mother 
particle: T = 0.01 Gyr (magenta line), 0.1 Gyr (green line), 1.0 Gyr (red line), and 10.0 Gyr (blue 
line), respectively. Black dotted lines with steeper slope are the fitting lines to the case with particles 
having a constant momentum, which are proportional to (1 + z)~ 3 / 2 , and the others are to the case 
with non-relativistic particle, which is proportional to (1 + z)~ 1//2 . See main text for details. 



constraints from a simple argument that the scale of any observed structure bounded by dark 
matter should be larger than the free-streaming scale of daughter particles. 

Free-streaming scale Ifs(z) at each redshift z can be estimated using the averaged 
velocity of daughter particles v(z) as 

l FS ( z )~v( z ) X -±— . (4.9) 

Here we define v(z) as 

rPth , . 

4irql > v(qT),z)fb > dqT) 



viz = 











Pth , m 
iTrqlfgUqv 



where 



<?d + mm - 2 



Here v(qj), z) is the magnitude of physical velocity of the daughter particle whose comoving 
momentum is (/d at redshift z. 

We calculate the free-streaming scale Zfs of daughter particles at each redshift z, which is 
shown in figure 5. We explain these curves in the matter dominated epoch (1 < z < 3000) as 
follows. First, when T^ 1 > H" 1 the daughter particles with a constant physical momentum 
Pth are kept being created and the averaged physical velocity of daughter particles becomes 
constant. Thus 

Const. , q/9 

FS ~17(yr ( } 

where we have used the fact that H(z) is proportional to a -3 / 2 = (1 + z) 3 / 2 in the matter 
dominated epoch. 
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Figure 6. Constraint on the lifetime of decaying dark matter from the free-streaming scale. The 
range painted over in gray corresponds to > H^ 1 (z = 3), at which the gravitational potential 
from the mother particles would be significant at z = 3, and hence it cannot be excluded. 

Secondly, when T _1 < H , i.e., the decay process has finished, the averaged velocity 
of daughter particles decays as v oc a -1 = (1 + z). Therefore 

/FS k { ~Wzf K (1 + Z) ~ 1/2 ■ (411) 

One can see that these dependences well describe the calculated curves in the figure. 

For the dark matter to form structures such as dark halos, the free-streaming scale Zps 
should be less than the size of the structures. We find in figure 5 that the free-streaming 
scale of daughter particles is sensitive to the lifetime and wiDi/ m M- Contrary to the 
case with massless daughter particles, the free-streaming scale of daughter particles becomes 
smaller if the lifetime of the mother particle becomes shorter. The reason is that the velocity 
of massive daughter particles decays faster in the earlier universe, because the expansion of 
the universe is faster. Hence the velocity of daughter particles decays in shorter timescale 
if they decay earlier, which leads to the smaller free streaming scale. As expected, the free- 
streaming scale becomes smaller if the mass ratio ttt-di/wm becomes smaller. Therefore any 
existence of large scale structure by dark matter can be used to constrain and m-Di/^-M 
through the free streaming scale. 

From the observations of Lyman a cloud at z < 3, the density fluctuations at about 
1 Mpc comoving scale have been found, for example, in SDSS [28, 36]. Therefore, when 
r _1 <C H~ x {z = 3) is satisfied, the range of mass ratio thdi/wim is excluded if Zps > 1 
Mpc (0.25 Mpc in physical scale) at z = 3. By taking this into account, we obtain the 
constraint on the lifetime of the mother particle as shown in figure 6. We find that the free 
streaming scale has a constraining power even for (1 — moi/V^Af) ^ 0.01, and the constraint 
is complementary to those obtained from the geometric distances to CMB and BAO. In the 
case T" 1 > H~ 1 (z), on the other hand, we can not use this method since most of the mother 
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particles, and therefore the gravitational potentials, still remain at the redshift z. A full 
treatment of cosmological density perturbations will be necessary in this case. Because our 
Lyman a constraint is given at z = 3, the region corresponding to > H~ l (z = 3) (grey 
region in figure 7) can not be excluded from the free streaming scale for now. 

4.4 Comparison with Peter et al. 



Figure 7. Comparison with the results from Peter et al. [7] on the life time and the kick velocity 
plane. Red, blue and yellow regions are the parameter space that we exclude in this paper. 

Before ending this section we compare our result with the constraints obtained by Peter 
et al. [9]. They considered a dark matter decay and describe it in terms of the velocity 
of daughter particles at their creation (they call it "kick velocity"). In a virialized dark 
halo, matter moves slower than its associated virial velocity iVir 

[37]. When > tVir, dark- 
matter halos will be disrupted by these particle decays. They performed several numerical 
simulations in order to study the detailed evolution of the total mass and density profile 
of galaxies composed of particles that undergo such velocity kicks as a function of the kick 
speed. As a result, is strictly restricted from the stability of the halos. We find that our 
constraints are comparable with their results in high kick velocity region, as shown in figure 
7. Note that our constraints are completely independent from theirs and we believe that our 
constraints are less uncertain in that they are free from the variety of galaxies. 

5 Conclusion 

In this paper we consider a decaying dark matter model in which the massive mother particle 
decays into two massive and massless daughter particles after cosmological recombination. 
We derive a complete set of Boltzmann equations to describe the evolution of the particles. 
We obtain constraints on the lifetime of the mother particle T _1 and the mass ratio w-di/^m 
with tod2 = from the Hubble parameter, CMB and BAO. The allowed T" 1 decreases 
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monotonically as mm/m-M increases. For the free streaming constraint, on the other hand, 
we find the opposite dependence. We find > 30 Gyr at the massless limit of daughter 
particles and mni > 0.97tom at the limit T^ 1 — > 0, from the distance to CMB (la). We also 
obtain constraints from the free-streaming of daughter particles from observations of Lyman 
a as (r-VlO- 2 Gyr) < ((1 - m m /m M )/W- 2 y 3/2 for T" 1 < H~ 1 (z = 3). However, to 
extend the free-streaming constraint to the range T" 1 > H^ 1 or to include the information 
from density perturbations such as CMB angular power spectrum, a complicated calculation 
of density perturbations is necessary. This will be an interesting subject and presented in a 
separate paper. 
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Appendix a Massless limit of daughter particle 

In the case that a mother particle decays into two massless particles, the energy densities of 
mother particle pm and daughter radiation p-Q satisfy the following equations 

Pm = —3HpM — TpM , 

and 

p D = -4,Hp D + Tp M ■ 
These differential equations can be solved to give 

PM = PM$a~ 3 exp(-ri) , 

and 

r PM0 f 



a 4 Jo 



where /9m0 is the expected energy density of dark matter without decay. In this calculation 
we obtained the following constraints, 

r^ 1 > 30Gyr , 

and 

T" 1 > 18Gyr , 

from the angular diameter distances to CMB and the position of BAO, respectively. These 
values agree with the results in the massless limit obtained from calculations in the main 
text, where the distribution functions of daughter particles are directly integrated instead 
using the above simple equations. 
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Next we show how one can obtain a limiting value of the mass ratio mni/mM in the limit 
r~ x —7- 0. In this limit our model corresponds to the A-CDM model such that Qm is replaced 
with — — f2 M where f2 M is the expected density parameter of dark matter without decay. 

We can deduce a constraint on the ratio mDi/fiM from a constraint on f^M from the angular 
diameter distances to CMB and BAO. In this way we obtain the following constraints, 

"mm/mm > 0.97, 

and 

"mm/mm > 0.23, 

from the angular diameter distances to CMB and BAO, respectively. These values are roughly 
in agreement with the values in figure 4, in which the distribution functions of daughter 
particles are directly integrated. 



Appendix b Present tempreture of CDM mother particle 

The temperature of mother particles at present time Tmo can be estimated by taking into 
account of the decoupling temperature of mother particles [38]. The temperature which 
corresponds to the mass of the mother particle mM = 1 TeV is so high that all species of 
standard particles such as eight gluons, W , Z°, three generations of quarks and leptons, 
and one complex Higgs doublet are relativistic at the decoupling of mother particles. We 
define g* as the total number of effectively massless degrees of freedom. In the epoch that the 
mother particles are in thermal equilibrium, we expect <?*( ca riy) = 106.75, while <7*( now ) = 3.36 
at present. 

The comoving number density of photons when the mother particles are in thermal 
equilibrium is 

9* (now) , v 

n 7 = n 7 o . (5.1) 

9* (early) 

Here the number density of photons at present n 7 o is 410 cm" 3 [36, 38] and therefore 

n 7 = 12.9 cm" 3 . 

On the other hand, the baryon number density can be derived from the baryon-photon ratio 
today 7] = n^/n^, which is estimated as rj = (6.19 ± 0.15) x 10 -10 (la C.L.) [39]. The ratio 
of the number densities between CDM mother particles and photons today, bm/^tOi can be 
written as 

= =n *-F> ' ( 5 - 2 ) 

n 7 o nb n 7 o 1 2b m M 

where fib is the cosmological density parameter of baryon. Here we have neglected a con- 
tribution from helium. By substituting rriM = 10 TeV, £2m = 0.222, and f^b = 0.0446 into 
eq.(5.2), we obtain 

^ = 2.89 x 1(T 12 . 

n 7 o 

The time of the thermal decoupling of mother particles can be estimated through the relation 
of the Boltzmann factor as, 

= exp — — 

n 7 V T d 
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Thus we obtain the temperature of mother particles as 

T M0 = a 2 d T d ~ 1.7 x 10~ 14 K = 1.4 x 1(T 18 eV , (5.3) 

where is the scale factor when the mother particles decouple from the thermal bath. We 
can write a d as 



where T 7 o = 2.725 K is the temperature of CMB 

Appendix c First order Boltzmann equations for daughter particle 

In this appendix we give a set of Boltzmann equations at first order which are necessary to 
compute the evolution of density perturbations associated with the mother and daugher par- 
ticles. The standard linear theory of density perturbations has been presented, for example, 
in the synchronous and the conformal Newtonian gauges [40]. Here we expand the linear 
perturbation theory by taking into account the decays of a SUSY CDM particle into two 
daughter particles [11]. 

By substituting eq.(2.10) and (2.11) into eq.(2.1) and (2.3), respectively, and comparing 
the equations order by order, we obtain the following equations at first th order 



1st order : * M + i-^-(h ■ h)^ M + n \ V ~ -{h + &h)(k ■ nf = . (5.4) 
aeu dmqu \ 2 



Here €m = y <?m + ° 2m M * s * ne energy of mother particles, and h and r\ are the metric 
perturbations in the synchronous gauge defined from the perturbed (space-space) metric hij 
as 

hij = J ykikjh(k,t) + (kikj — —6ij)6ri(k,t) \ exp(—ik ■ x)d 3 k , 

where k = kk is a wave number vector. Hereafter, we omit the arguments of h(k,t) and 
r](k,t) for simplicity. For daughter particles, we obtain 

din f / 1 

1st order : i> Dj + i^~(k ■ n)^ Dj - + ( fj - - (h + 6r))(fe • n" 2 

1 



aej)j dln^D \ 2 



AO) 

J-Dj 



r(gD,9M)/M ) (gM)(^M(gM) - ^Dj{qD))d 3 qM 

Note that on the RHS of the above equation, while the first term directly comes from the 
collision term, the second term results from the change of the background evolution of the 
distribution function, i.e., /P.. 

Following [40], we consider the Legendre expansion of and \£d into *&mi and ^>bi, 
respectively. Here the coefficients VPmz and ^dz are defined as, respectively, 

+00 

tf M (fc, n, q, t) = J^(-i)'(2/ + l)*Ml(fc, q, t)Pi(k • h) , (5.6) 
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and 

+00 

^ Dj (k,h,q,t) =Y,(-i) l (^ + ^B ] l(k,q,t)P l (k ■ n) , (5.7) 
1=0 

where Pi(k ■ h) are a series of Legendre polynomials. The factor {—if {21 + 1) is chosen to 
simplify the expansion of plane wave. Note that equation eq.(5.4) is the same as the case 
of massive particles such as neutrinos without decay [40]. Therefore we obtain the recursion 
formula for the mother particles without any explicit V dependence as [40], 

9M^, T , , 1 L dln fM 



^mo = *Mi + 7;h , 

aeu o dmqu 

*M1 = |^(#M0 - 2^M2) , 
oaeM 



*M2 = i (2*M1 - 3*M3) - \-^ h + TV 



(0) 



ClM k /n,T, a.T, \ ( 1 I , 2 A dln fM 



5aeM \15 5 / dlnqM 

*M(n) = / m ^ («*M(n-l) ~ (« + l)*M(n+l)) , ( n > 3 )- 

v ; {In + l)aau 

Hereafter, we set \&M(n) = f° r n > 2, because we assume that the mother particles are 
CDM. For daughter particles, we obtain 



A _ 1»k l ^ln/W 4r ^ 0) (^) ^ , » 

ae D j 6 dlng D #M <?d f£>(q D ) 



3ae Di ^M^d/d («d) 

- J^T^^in^ ~ (» + - £i^* Diw(flD ) (n > 3). 

In the limit T — ?■ 0, the above equations reduce to those of massive neutrinos in the standard 
perturbation theory [40]. Again, the terms which are related with the decay process comes 
both from the collision term and the change of the background evolution. For the moment 
equations of daugher particles with n > 2, the decay terms only come from the change of 
the background evolution, because the mother CDM particles do not have the corresponding 
moments, i.e. ^M(n) = for n > 2. 
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